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Abstract. We consider a three dimensional Wigner crystal of electrons lying in a host ionic dielectric. 
Owing to their interaction with the lattice polarization, each localized electron forms a polaron. We study 
the collective excitations of such a polaronic Wigner crystal at zero temperature, taking into account 
the quantum fluctuations of the polarization within the Feynman harmonic approximation. We show that, 
contrary to the ordinary electron crystal, the system undergoes a polarization catastrophe when the density 
is increased. An optical signature of this instability is derived, whose trend agrees with the experiments 
carried out in Nd-based cuprates. 

PACS. PACS-key describing text of that key - PACS-key describing text of that key 



1 Introduction 

Let us consider a Wigner crystal (WC) of electrons — a 
state where the carriers form an ordered array in order 
to minimize the repulsive Coulomb energy — and put it 
into a polar dielectric — a material which can respond to 
the electron motion by displacing the positive and nega- 
tive ions. If the density of charges added to the dielectric 
is low, each carrier independently forms a polaron, i.e. it 
is surrounded by a polarization cloud due to the interac- 
tion with the longitudinal optical phonon modes of the 
host medium. In the opposite high density limit, dynami- 

\ cal screening will prevent both the localization of charges 
due to the Coulomb repulsion and the formation of pola- 
ronic bound states, thus leading to a weakly interacting 
liquid. The phase transition occurring upon varying the 
density between these two extreme limits has been anal- 

\ ysed in ref. [[j] starting from the crystallized state. The 
resulting picture allows to distinguish between two qual- 
itatively different situations. While in the weak and in- 
termediate electron-phonon coupling regimes, the polaron 
Wigner crystal (PWC) can melt towards a polaron liquid, 
this is prevented when the polar coupling is strong: in that 
case the polarization cloud is virtually frozen and polarons 
tend to dissociate as the system becomes metallic. 

The scope of this work is to go beyond the mean-field 
approach presented in ref. |l]] , with particular attention to 
the collective modes arising from the long-range interac- 
tions between electrons in the crystallized phase. The main 
result of our study is that in the strong electron-phonon 
coupling regime, where polarons have a well-defined in- 
ternal structure, the dipolar electron-electron interactions 
can lead to very anomalous dielectric properties, and even- 
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tually to an instability at high enough density, character- 
ized by a softening of a long wavelength transverse col- 
lective mode. An equivalent result was already obtained 
in previous works assuming a static polarization poten- 
tial [Q], i.e. neglecting the quantum fluctuations of the 
polarization. Similar indications have also been given by 
other authors || starting from a classical liquid state. All 
these cases are practical realizations of the very general 
phenomenon of a "dielectric catastrophe" 0), occurring 
when local dipoles oscillating with a given restoring force 
(the electrons inside the polarization potential-wells) in- 
teract through the long-range Coulomb forces, which tend 
to soften the transverse restoring potential. Such instabil- 
ity, which can be traced back to the local Lorentz field, 
has an experimental signature: the frequency of the po- 
laron peak in the optical conductivity decreases with dop- 
ing, a trend that has indeed been observed in underdoped 
Nd2- x Ce a: Cu04_j / (NCCO) below the insulator-to-metal 
transition ||. 

The paper is organized as follows. The main features 
of the model, which was thoroughly introduced in refer- 
ence H, are summarized in section ^. In section we 
calculate the phonon dispersion of the PWC on the ba- 
sis of Feynman's harmonic trial model for the polarons. 
In section [|, we evaluate the effect of the vibrations of 
the PWC by calcultating the dielectric constant e (k, oj) 
as a function of the density. We show that it is negative 
on a large region of k and u> when the system approaches 
the instability. In conclusion, we propose a scenario where 
free carriers liberated by the dissociation of polarons could 
coexist with the residual localized polarons in a "mixed" 
phase beyond the instability, and speculate about a novel 
pairing mechanism where the collective vibrations of the 
PWC or, equivalently, the "overscreening" related to the 
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negativity of the dielectric constant s(k,u>), acts as a glue 
for the Cooper pairs. 



2 Model and approximations 

The model for the crystallized state has been thoroughly 
defined in ref. equations (f)-(5). Electrons with band 
mass to* at a density n = (47r/3) 1 R~ 3 interact directly 
through the long-range Coulomb forces. These are im- 
mersed in a polar background characterized by the longi- 
tudinal phonon frequency ujlo and the effective dielectric 
constant e. The polar medium responds to the electron 
motion through a Frohlich type interaction characterized 

1/2 

by a coupling constant a — (m* /2h ujlo) s 2 /e, which 
can lead to polaron formation. Once the degrees of free- 
dom of the polarization field are integrated out, the re- 
sulting action for the electrons can be expanded assuming 
small displacements relative to the equilibrium posi- 
tions. To quadratic order, this reads: 



S({ Ui })=5> + i5> ; 



with 
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where the indices a, 7 = (x, y, z) denote the cartesian co- 
ordinates. The dipolar matrix elements are defined as: 



A ai = 

13 



5 a ^ R^j 3i?j j RJj 



(4) 



We have introduced the phonon propagator G w (/?, r — a) = 
(n+\)e-^ T - a \ +ne^ T - a \ together with n = (e^-1)" 1 . 
The frequency uj\y is related to the plasma frequency by 
the relation w 2 ^ = w 2 /3 = e 2 /m*R 3 . (5 is the inverse tem- 
perature, and the units are such that h = ks = 1. 

Expressions P)-(|^) constitute the basic model for the 
polaron Wigner crystal (PWC). It is based essentially on 
the following approximations: (i) anharmonic terms in the 
expansion in {u^} are not considered; (ii) exchange ef- 
fects between electrons are neglected, and so is the possi- 
ble (weak) magnetic ordering of the crystal. For the main 
purposes of our studies, i.e. the crystal stability, all these 
approximations are well justified in the range of densities 
that we consider, as was discussed in detail in ref.[0. 



The local part of the electron action has been exten- 
sively studied in ref. Q , by isolating neutral entities con- 
taining only one electron — the so-called Wigner spheres 
- which are non-interacting at mean-field level. In this 
approach, all the correlation effects are carried by a local 
harmonic potential [the second term in square brackets in 
eq. @], which merely adds to the polaron problem (the 
additional constant term is the electrostatic energy com- 
ing from the jellium approximation). The last term in eq. 
(||) is the self-induced interaction mediated by the po- 
larization, which is responsible for polaron formation. It 
is not solvable, but it can be treated successfully within 
the Feynman variational approach || by introducing a 
quadratic trial action So of the form 



. 2 m* Jly 
u dr — 

2 e s 



L 2 dT 



^ I I \u{t)-u(s)\ 2 G w {/3,T-<T)dTda, (5) 



(we have dropped the site index, since at this stage all the 
sites are equivalent). It can be seen that expression (0) 
results from the path-integration over X of the following 
Lagrangian: 

to* . 9 M ■ n K , _,.2 m* ujIv , , . 
Lo = ~u 2 + -X 2 - y (u-X) 2 --^V- (6) 

where w 2 = K/M. In eq. (|6|), K and w are variational 
parameters to be adjusted in order to minimize the free 
energy. This method of approximation, not only yields the 
best known analytical results for the polaron energy, but 
it also gives a very good insight into the physics of the 
polaron problem: when the variational parameters K and 
w are properly adjusted, the model @ is a good schema- 
tization of the true polaron, especially at low tempera- 
tures (roughly at T < lo^o, where the polarization cloud 
is sufficiently "rigid"). Basically, the actual polarization 
of the polaron is replaced by an auxiliary particle with 
coordinate X and mass M , which reduces the many-body 
polaron problem to a two-body problem. 

Let us now consider the non-local part (||) for i =/= j. 
These terms, which represent dipole-dipole interactions 
between different electrons, become important beyond the 
mean-field approximation described above, and they are 
responsible for the dipersion of the vibrating collective 
modes. The first term in eq. (||) is instantaneous, and it 
comes out from the direct Coulomb electron-electron re- 
pulsion. The second term is a retarded dipole-dipole inter- 
action which, as the polaron term in (S), is mediated by 
the polarization. Its physical meaning is that an electron 
i, when it moves, feels the polarization field created by 
electron j. It is important to emphasize that, unlike the 
polaron term in (||) , this retarded dipole-dipole interaction 
vanishes with uj^o (he. a — > 00): when the polarization is 
completely frozen, only the instantaneous dipole-dipole in- 
teraction remains at work. This particular limit, which is 
very instructive, has been already studied in ref. Q. Here, 
we consider the general case ujlo 0; where the retarded 
dipole-dipole interactions cannot be neglected. 
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2.1 Effective quadratic model 

In the same spirit as the Feynman approach for the single 
polaron, we introduce an effective Lagrangian: 



n* , 9 M ■ 2 K , ,^ . 2 r, 
Y Ui + Y Xi -2^-^—2 e 



L eff = E 

i 



IV 2 



(7) 



where the Ay and By are 3x3 matrices in coordinate 
space. This model can be thought of as a generalization 
of the Lorentz lattice of dipoles, which includes retarded 
dipolar interactions through the auxiliary particles Xj. 

The effective action for the electrons is obtained after 
path- integration on the coordinates X,-: 



s eff - 



J2 [ dr[vZ(r) - (K/m*+^ w /e s )uUr)} 

^E / u »( T )4 u j( r ) dr 

dTdaG w ((3,T — a) x (8) 
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At this stage, the most rigorous way to proceed would 
be to self-consistently solve the full Feynman variational 
problem given by: 



E = min[_En 

K,w 



lira < S — S e f / > 



(9) 



where Eq is the zero-point energy of all the eigenmodes of 
(^) However, we can get more insight into the problem 
by making use of two additional approximations: 

(i) We restrict to a single variational parameter, by 
putting w — ujlo in the expression (JsJ) . As was pointed 
out in ref. [[TJ , this has very little effect on the results when 
dealing with moderate to large e-ph couplings (a >~ 5). 
From now on, we will mainly focus on this case, which is 
to our opinion the most interesting one. 

( ii) For the variational parameter K, we use the results 
obtained in the framework of the mean-field approxima- 
tion (see ref. Q). The underlying idea is that of a pertur- 
bation expansion around the mean-field result, so that the 
first correction is calculated by taking the zeroth order as 
a good starting point. At low enough density, the dipo- 
lar terms are indeed a small perturbation, since they ex- 
plicitely involve the fluctuations of the electrons around 
their equilibrium positions, which are small in this limit. 
More quantitatively, in the case of the ordinary WC, it was 
shown that the inclusion of dipolar interactions changes 
the total energy by no more than 5-10 percent [pi. 



Of course, the previous argument no longer applies 
when approaching the critical density, where one can ex- 
pect some quantitative deviations from the predicted phonon 
spectrum. However, our results should be at least qualita- 
tively correct in all the insulating phase. 

Within this framework, there is no need to perform 
new path-integrations: the matrices Ay and By are cho- 
sen in order that all the dipole-dipole terms of S — S e f / 
in equation (||), vanish term by term. This leads to the 
following set of equations: 



A, 



-2KB U + J2 BliBu = 



-Ay ; i + j (10) 

1 

(11) 



-2KB, + B * B u = K& -h ; * ^ 3- (12) 



The first equation straightforwardly defines the matrix Ay 
of the instantaneous dipole-dipole interactions, and the set 
of equations for the retarded terms By will be solved in 
the next section. As a result, we are able to calculate the 
phonon spectrum of the effective Lagrangian (|7|), which in 
the same manner as (^|) for the single polaron problem, 
mimics the collective excitations of the PWC. 



3 Phonon Spectrum and Instability 

In this section, we derive the phonon spectrum of the La- 
grangian (0), as was done for the usual Wigner Crystal 
and Lorentz lattice of dipoles in references p[|l0|,[8| . 

First, the WC of electrons is recovered from eq. (0) by 
putting the e-ph coupling to zero, with e s — * Eoo- All the 
terms involving {X^} thus vanish, and we obtain [pjj|: 



V ^ TTl* . n ™* Ww 2 1 a 

L w = l^— u i~—^r Lu i~o E u i A n u i 



(13) 



with the matrices Ay for i =/= j given by equation (|10|). 
The potential term of ( |l3| ) can be included in the defi- 
nition of the A's by putting An = (m*uj^ v /£ 00 )I, with / 
the identity matrix. We can now diagonalize the coupled 
equations of electron motion 



by introducing eigenmodes of the form: 



/,n i(k-Rjj — 

U;(t) = UkA^kAe 







(14) 



(15) 



where ekA is the polarization vector of the mode (k, A), k 
being a wavevector in the Brillouin zone, and A the index 
of the (three) phonon branches. If we define the Fourier 
transform as 



Ax, 



N ^ 



(16) 
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then the eigenfrequencies fulfill: 

Ar3 *w 2 (k,A) 
ek\NA k e ky = to* 



5\\> 



(17) 



For a given Bravais lattice, the matrix elements of A k are 
given by [f|: 

A£» = k"k v V k + V k+K {k» + K v ){k» + K") 

-Y^K^K v Vk. (18) 

where 14 = 47re 2 /£; 2 £ 00 is the Fourier transform of the 
Coulomb potential, with Vk=o set to zero for charge neu- 
trality, and K a reciprocal lattice vector. From (fL7h and 
( |i~8| ) it is easy to prove the following result: 



w (k,A) 2 



A 



(19) 



which is known as the Kohn sum rule H[LC]]. The phonon 
spectrum of the WC contains two acoustical transverse 
modes (t) and one optical longitudinal mode (£), which 
satisfies uj(0,£) — Wp/^/ioo. 

Let us now calculate the Fourier transform £?k of the 
matrices Hy, defined in the same way as A k in eq. (|l6|). 
Multiplying eq. (0) by cxp(zk-Ry ■), summing over j, then 
making use of (|1 1|) , one gets: 



2KNB V + NB k NB u 



K—NAk 

£ 



I (20) 



The above equation can be solved in the same basis {ek.A} 
which diagonalizes ^lk and A^. Therefore, one obtains 
ek,A-^V-B k e k ,A/ = 5(k, A)<5 a ,a/ where: 



3.1 Phonon spectrum 



1 to* w 2 (k, A) — 



e M 



J LO 



(21) 



We can now determine the whole phonon spectrum of the 
PWC by diagonalizing the Lagrangian ([?]). Let us intro- 
duce the Fourier transforms 



, WkA^kAe 



-ik-Ri 



x ^7i|? XkAekAe " lk ' Rl - 



(22) 



kA 

then the equations of motion become 

.2 , ,2 



"kA 

ATkA 



e e 
[K - g(k A)] 
M 



w 2 (kA) K 

TO* 

K 



, [K-g(KX)] v 
MkA H : -*kA 



-UkA 



M' 



kA- 



(23) 



fi(k,X) 



(Opal 
tOro 



(£ s Mp/m*)" 2 



continuum 
(free electrons) 



int. modes 
(Lorentz lattice) 



ext. modes 
(Wigner lattice) 



Fig. 1. The excitation spectrum of the PWC. 



For each wavevector k and polarization A, the solution of 
the secular equation gives the eigenfrequencies /2(k, A) of 
the PWC in terms of the eigenfrequencies w(k, A) (formula 
(18)) of the WC of electrons: 



(24) 



±t 



, 2 (k,A) 



" Mo' 



, 2 (k,A) 



where we have defined 



H^ + iKo-^ (25) 



J pol 



In the following, we shall denote by ]?j n t(k, A) [{] ext (k, A)] 
the solution with the + [— ] sign. By looking at the equa- 
tions of motion (p3|), one can verify that in the limit of 
low k, these frequencies are associated respectively with 
the out-of-phase and in-phase vibrations of the electrons 
relative to the polarization. Interestingly enough, within 
our approximation scheme, the polaron effect now enters 
in the frequency spectrum through the parameter u} po i 
alone. Moreover, the parameter oj po i is identified as the 
frequency of the transverse mode i?j n t(k, t) at k = 0. 



It is instructive to work out the solutions of eq. 
in the low density limit, where u>w is small compared to 
the internal frequency of an isolated polaron. In this case 
we obtain 



^L(XA) 



^L(k,A) ~ uj 2 pol + 



u 2 (K X)/e s 
Mp/m* 

w 2 (k,A) 



(26) 
(27) 



Those frequencies are sketched in figure |]. At low density, 
there is a clear separation of energy scales between: 
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i) the low frequency "external" modes with f2 ext <C 
UiO) which correspond to a Wigner crystal of particles 
with mass Mp instead of m* . At such low frequencies, the 
polarons behave as rigid particles and the internal elec- 
tronic motion can be neglected. Furthermore, all the po- 
larization cloud participates in the screening of the dipolar 
interactions, so that the eigenfrequencies are reduced by 
a factor e a . 

ii) the high frequency "internal" modes with flint 3> 
UiO, which correspond to a Lorentz lattice of dipoles with 
characteristic frequency u> po i, which was extensively stud- 
ied in ref. 0. In this case, the electron vibrations are very 
fast and the appropriate dielectric constant is £oo. 

Such separation of energy scales is strictly valid only at 
very low densities. At densities such that lo po i approaches 
lolOi i- e - close to the phonon instability (see below), the 
full expression (24) must be used to determine the eigen- 
frequencies. 

The treatment of this section provides us with an es- 
timate of the low-lying excitations of the PWC. At this 
stage, it is useful to comment on a few points concern- 
ing our approximations. First of all, above some threshold 
frequency higher than co po i , the calculated eigenmodes will 
eventually merge into a continuum of states, correspond- 
ing to free electrons ionized from both the polarization 
potential and the periodic potential of the Wigner crystal 
itself. In addition, more bound states can appear below 
such threshold, if there are metastable interstitial states 
with long enough lifetimes. As far as the Feynman approx- 
imation is concerned, we know from the single polaron 
results Q| that the interpretation of the eigenfrequency 
of the quadratic model as representative of the excitation 
spectrum of the real system is qualitatively correct in the 
strong coupling limit, where the optical Franck-Condon 
transition is very sharp (indeed, the polarizability of a 
system with a definite quantum transition has exactly the 
same form as for a classical harmonic oscillator). At in- 
termediate values of a, however, the transitions inside the 
polarization potential- well are considerably broadened. In 
this case, the calculated eigenfrequencies should still be of 
significance, provided that we associate to them a finite 
linewidth r, related to the so-called Huang-Rhys factor. 



3.2 Instability of the polaron system 

Contrary to the ordinary WC of electrons, the PWC can 
undergo a phonon instability when the density reaches a 
critical value. This can be most clearly seen by writing 
the action (JsJ) as a sum over the quadratic actions of the 
independent eigenmodes UkA, Seff — SkA ^kA, with 

SkA = -^ / |^ k A(T)| 2 rfT~ m *^ (k ' A) f \u kX (T)\\M) 



K(k, X)lo lo 



\uk.x(r) - u k \(o-)\ 2 Gf3^ LO dr da 



term (the retarded interaction of Uk\ with itself) is now 

w 2 (k, A) - li>w 



K(k, A) =K 



1 m* 



to 



LO 



Both oj 2 (k, A) and are proportional to the density, so 
that for n — > we recover K(k, A) = K. However, upon 
increasing the density, the restoring force K (k, A) can be- 
come negative, and the system is unstable. The first insta- 
bility occurs due to the long wavelength transverse modes 
for which uj(k = 0,t) — 0. It is thus the long wavelength 
transverse vibrations of the electrons relative to their po- 
larization cloud which cause the instability, a phenomenon 
which is reminiscent of the phcnomcnological dissociation 
criterion introduced in reference Q. 

In terms of the parameter u) po i, we see that K{k = 0, t) 
vanishes (i.e. the system is on the verge of the instability) 
when 

u P oi = u LO (29) 

This instability condition is the correct generalization of 
the condition ui po i = of ref. |2| to the case of a fi- 
nite phonon frequency. Let us stress once again that the 
phonon instability occurs at k = 0, where the properties 
of the system are independent of the details of the lattice 
ordering, provided that the system remains isotropic. As 
a consequence, the instability criterion ( p9| ) also applies 
to disordered insulating polaron phases, which can be re- 
alized in real solids where the doping ions act as impurity 
potentials, thus driving the particles away from their con- 
figurations on a Bravais lattice to form a Wigner glass. 
The above criterion should also hold in the case of a po- 
laron liquid, provided that the motion of the polarons is 
slow compared to the vibrations of the electrons inside 
their potential wells pi[. 



a 


(melt ) 
'8 


1 s 


( inst ) 





960 






3 


537 




19 


7 


52 




18 


10 




32 


16 


100 




31 


15 



We see that each of the SkA has the same form as the 
mean- field action (||), but the coefficient of the polaron 



Table 1. The density parameter ri* n3t ' for the phonon in- 
stability, compared with the critical r s for crystal melting 
and polaron dissociation in the mean field approximation (see 
ref. M). Parameters are Eco — 5, e s = 30 and m* = 2m e . 
The critical densities can be obtained from the relation n = 
1.6 • W 24 /r'icm~ 3 . 



In Tabic I we have reported the instability parameter 
r (mst) ag deduced from criterion (^9|), compared to the val- 
ues for polaron dissociation and crystal melting obtained 
in reference Q. In the weak coupling limit, the retarded 
interactions can be neglected and the PWC tends to a 
WC of electrons, which has no phonon instability [fi"o|| . 

Correspondingly, ri mst ^ vanishes (the critical density di- 
verges) as a — » 0. In the opposite strong coupling limit, 
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3 102" 4 102" 5 10-" 6 10» 7 10 2 ° 
n (cm- 3 ) 

Fig. 2. Softening of the collective frequency ui vo i with density. 
The instability occurs when uj po i = ujlo- Parameters are a — 
7, Sao — 5, e s — 30 and m* = 2m e . 



u>j,q), lolo and loto being respectively the longitudinal 
and transverse optical phonon frequency, the force acting 
on an electron at site i can be written as 



P i (i) = -eE eBt (R 1> t) 



4-7rk 



e P0 i(k-r-o)t) 



k 2 



- ph 



(31) 



Introducing this driving force into the equations of mo- 
tion ( p3[ ) and solving for the induced displacement u k A we 
obtain 



iiirepo (k • e kA ) 2 



LO 



e ph (u>) k 2 [nf nt (k, A) - w a] [^ xt (k, A) - u 2 



The longitudinal dielectric constant can be obtained by 
comparing the total charge associated to the displace- 
ment MkA a- n d the perturbing charge distribution, namely 
e L (k,u) = Pext/ptot (for more details, see ref. @,||). The 
result for the PWC reads 



where ujlo — * 0, we recover the results of ref. [0. Gen- 
erally speaking, in all the range a > a* , the results for 
r (mst) a g ree w jtlvi n a factor of two with the critical r s 
for polaron dissociation, which confirms that there is an 
intimate connection between the instability of the long 
wavelength transverse phonons and the mechanism of po- 
laron dissociation studied in previous references. On the 
contrary, for a < a*, the crystal melts due to the quan- 
tum fluctuations of the polarons before the occurrence of 
the phonon instability, and the latter is never realized in 
practice. 

Figure || shows the evolution of ui po i vs the density in 
an intermediate situation, a — 7. The behaviour is very 
similar to the case lulo — presented in ref. j^j, where 
it was shown that the frequency of the collective mode is 
gradually reduced due to the effects of the local Lorentz 
field (see also ref. UM). 



4 Longitudinal Dielectric Constant and 
Optical Response. 

In this section, we derive the screening properties and the 
absorption spectrum of the PWC in the harmonic approx- 
imation. The determination of the dielectric constant of 
the model (Q) follows the lines given in ref. (pJ), and the 
optical conductivity is deduced from the k = limit of 
this quantity. 



4.1 Dielectric Constant 

The dielectric constant is defined as the response to a 
perturbing external charge, that we take of the form 



1 



1 



Pext(r,t) = poe^-^ 



Taking into account the screening of the host polar lattice, 
which we schematize as e p h{oj) = £oq{uj 2 — u>\q)/{uj 2 — 



Ei(k,w) e p h(v) 



E 



(k ■ e kA ) 2 



(32) 



[fl? nt (k,A)- W 2] [^(k,A)-, 



This expression can be simplified in the following limits: 



4.1.1 long wavelengths (k = 0) 

Since in this limit k • e k t = 0, the transverse modes do 
not enter in the sum over A, and the dielectric constant 
reduces to 



£i(0,w) = s p h(uj) 



(33) 

This result is sketched in figure @. In this limit, the trans- 



100 




(30) Fig. 3. Dielectric function of the PWC at k = 0, with the 



same parameter set as in figure Q at n = 5 ■ 10 cm 



From 



the Lyddane-Sachs- Teller relation, the transverse phonon fre- 
quency is loto = 0.41u;lo. 
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verse and longitudinal dielectric constants become equal, 
and the zeros and poles of (|3^) correspond respectively to 
the longitudinal and transverse collective modes. 

The pole at u = corresponds to a transverse slid- 
ing mode of the PWC, which is equivalent to the sliding 
mode of a CDW 0] . This mode is expected to be current- 
carrying when the PWC is not pinned, as is the case in 
our jellium model, leading to an infinite d.c. conductivity. 
In real materials, however, this excitation can be shifted 
to higher frequencies due to the pinning by impurities and 
commensurability with the host lattice. As a consequence, 
the infinite d.c. conductivity is desrtoyed and the infinite 
(negative) dielectric constant at u = is replaced by 
a large (positive) value below some threshold frequency 
Api n characteristic of the pinning process (this effect can 
be introduced in our model in a phenomenological way by 
taking the frequency Q ex t (k, oj) of the transverse mode to 
be equal to A p i n at low k). 



4.1.2 static limit (oj = 0) 



In this limit, eq. 



becomes 



Ei(k,0) 



k 2 ^ 



(k ■ e kA ) 2 
^ 2 (k,A) 



(34) 



By comparing this expression with the result of ref. |10| , 
we see that the static response of the PWC in the jellium 
model is identical to the ordinary WC, with an enhance- 
ment e s due to the host phonon screening. 



4.1.3 sign of the dielectric constant 

By taking a simple model for the phonon dispersion of 
the WC of electrons, we can study the sign of the longitu- 
dinal dielectric constant of the PWC in the (k, oj) space. 
The model consists of two transverse acoustical branches 
ct>(k, t) = ak and one longitudinal optical branch w(k, I) = 

in order to fulfill the Kohn sum rule. 



oj 2 - 2a 2 



An 



estimate of the parameter a (the velocity of sound) can 
be obtained by arguing that the value of the transverse 
frequency at the end of the Brillouin zone (ko = n/R s ), 
must be equal to some fraction of the plasma frequency 
(we shall take a — 0.1oj p /ko). Substituting these expres- 
sions into eq. ( f24|) we get the eigenfrequencies of the PWC, 
that we use in eq. (p3) to determine the dielectric con- 
stant as a function of the density. For simplicity, we re- 
strict k to vary along one of the principal axes of the 
crystal, where by symmetry the transverse polarisation 
vectors obey e k t ■ k = 0. 

As can be seen in figure g, the regions where El < 
expand with increasing density. When one approaches the 
critical density, the dielectric constant is negative in all 
the region below the frequency Qi n t(k,£) ~ oj p / ^/£oo- 




o 
j 
3 
3 



0.4 0.6 

k/(7i/R s ) 

Fig. 4. Regions in the (k,ui) 
stant is negative, which are 
the poles of el- Left panel: 



0.2 



0.4 0.6 0.8 
k/(7C/RJ 



plane where the dielectric con- 
located between the zeros and 
n — 5 ■ 10 19 cm~ 3 ; right panel: 



n = 1.7 ■ 10 cm . Parameters are the same as in figure 1. 
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Optical conductivity of the PWC at various densi- 



Fig. 5. 

ties: n — 0.5, 1.7, 2.3 and 2.4 • 10^ u cm J , with the same pa- 
rameter set as in figure |^ (the instability occurs at n' lnsi ' = 
2.7 ■ 10 20 cm -3 ). We have introduced the phenomenological 
linewidths r po i — ujlo, i~b — r p h = Q.Iuolo- Note the shift 
of the polaron peak, and the transfer of spectral weight taking 
place between the feature at ui vo i and the zero frequency peak. 
When the polaron lattice is pinned by impurities or by com- 
mensurability effects, the u) = sliding mode moves to a finite 
frequency to = A p i n . 



4.2 Optical Conductivity 

The real part of the optical conductivity is related to the 
imaginary part of the transverse dielectric constant by 



UJ 



o(w) = -7— I'm e(k 
4tt 



0,w) 



(35) 



We remark that expression (g^) is perfectly real, and there- 
fore it can only lead to delta-like absorption peaks located 
at its poles uj = and u> = u> po i . As was stated in the pre- 
vious section, at finite values of the e-ph coupling, the col- 
lective excitations have an intrinsic lifetime which cannot 
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be obtained within the quadratic model (|7|) . We thus intro- 
duce phcnomcnological lincwidths by putting lo — ► LO + iT. 
Substituting the values of L?.;„ t (0, A) and Q ex t(Q, A) given 
by eq. (Jm)), we obtain 



(36) 



where the first term is the response from the phonons of 
the host medium, the second term comes from the collec- 
tive transverse excitation at lo po i , and the last term is the 
contribution from the sliding mode centered at lo — 0. As 
was stated in the previous section, if the PWC is pinned 
due to disorder or lattice commensurability, a gap opens 
in the optical conductivity, and the zero-frequency peak 
is shifted to u> = A P i n . Therefore, the optical response of 
the polaron Wigner crystal is characterized by the emer- 
gence of two collective-mode peaks, a polaron peak lo- 
cated above the host phonon frequency and a sliding (or 
pinning) peak around or below the phonon frequency. 

Assuming that the above contributions are well sepa- 
rated from one another (i.e. that the peaks do not overlap), 
these take the form of Lorentz peaks: 



tTph(w) = ; LO, 



4-7T 



TO 



(u>-UTo) 2 + I? h /i 



r/4 



ct (lo) 



Upot) ^(to-Lo po i) 2 + r 2 /4 
'uj 2 lo \ lo 2 r 



(37) 



(38) 



(39) 



The total spectral weight associated with the absorption of 
the PWC is lo 2 /8, and it obeys the conductivity sum rule. 
At low density, if a > a*, one typically has lo po i 3> lolo 
and almost all the spectral weight is carried by the po- 
laron peak in a po i{io), as can be seen from eq. (|38|). As the 
density is increased, some spectral weight is transferred 
to the sliding mode at lo = (lo — A p i„ in the pinned 
case), until the polaron peak eventually disappears at the 
instability point. 



5 Conclusion 

In this paper, we have analyzed the melting of a polaron 
Wigner crystal in the strong electron-phonon coupling limit, 
taking into account the dipole-dipole interactions between 
localized electrons, i.e. going beyond the mean field ap- 
propach developed in ref. M. We have shown that the 
system undergoes a phonon instability due to the effect 
of the Lorentz local field, which confirms the polaron dis- 
sociation scenario introduced previously. The instability 
has an optical signature: the polaron peak in the optical 
conductivity shifts towards low frequencies as the doping 
density is increased. This red-shift should be accompanied 
by the rise of a collective peak in the far infrared, when- 
ever the polaron Wigner crystal is pinned by disorder or 
by commensurability effects. This complex behaviour has 
been experimentally observed in NCCO ||, which sug- 
gests that both the polarons and the unscreened long-range 



Coulomb interactions should play an important role in the 
high Tc compounds. We suggest that the same behaviour 
could be observed in other superconducting materials. 

We have calculated the dielectric constant of the PWC, 
and shown that it is negative on a large region of (k,Lo). 
This indicates the possibility of reaching a superconduct- 
ing phase, if carriers of different nature — localized po- 
larons and free electrons — can coexist close to the in- 
stability density. In such a regime, the itinerant electrons 
would be paired through the high frequency internal vi- 
brations of the polaron lattice, which would act as the 
phonons in the ordinary BCS theory The large collective 
frequencies of the polaron crystal would naturally lead to a 
high superconducting transition temperature. Eventually, 
at higher doping levels, the polaron clouds are completely 
screened by the free electrons, so that the mixed phase 
(and the superconducting phase) should disappear in fa- 
vor of a more conventional metal. Jll| 

It should be kept in mind that our results are based 
on the continuous Frohlich model in the low density regime 
and at strong electron-phonon coupling, which misses many 
essential physical features of the cuprates. Indeed, we are 
not considering the spin degrees of freedom, although these 
are certainly responsible for most of the anomalous prop- 
erties of the normal state. We are not taking into account 
the formation of stripes, although the practical realiza- 
tion of the mixed phase conjectured above could well in- 
volve such form of microscopic phase segregation. On the 
other hand, the inclusion of other ingredients, such as the 
anisotropy of electronic bands in the cuprates, or the dis- 
cretization of the lattice, would increase the tendency to 
localization and are therefore expected to push the system 
towards the present "polarization catastrophe" scenario. 

The authors would like to thank P. Calvani, B.K. Cha- 
kraverty, S. Ciuchi, D. Feinberg, D. Mayou and J. Ran- 
ninger for many stimulating discussions. 
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